System-reservoir dynamics of quantum and classical correlations 
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We examine the system-reservoir dynamics of classical and quantum correlations in the decoher- 
ence phenomenon within a two-qubit composite system interacting with two independent environ- 
ments. The most common noise channels (amplitude damping, phase damping, bit flip, bit-phase 
flip, and phase flip) are analyzed. By analytical and numerical analyses we find that, contrary to 
what is usually stated in the literature, decoherence may occur without entanglement between the 
system and the environment. We also show that, in some cases, the bipartite quantum correlation 
initially present in the system is completely evaporated and not transferred to the environments. 
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I. INTRODUCTION 

Until recently the quantum aspects of correlation were 
attributed to inseparable quantum states [1], i.e., all non- 
classical correlations in a composite quantum state was 
regarded as entanglement. However, the discovery that 
mixed separable (unentangled) states can also have non- 
classical correlation [2, 3] and that the use of such states 
can improve performance in some computational tasks 
(compared to classical computing) [4, 5] has opened a 
new perspective on the study and comprehension of such 
correlations. The distinction between quantum and clas- 
sical aspects of correlation in a composite quantum state 
is an important issue in quantum information theory 
(QIT). It is largely accepted that quantum mutual in- 
formation is a measure of the total correlation contained 
in a bipartite quantum state [6, 7], but an outstanding 
question is how to distinguish between the quantum and 
the classical aspects of the total correlation. In view of 
the distinct nature of the correlations (quantum and clas- 
sical), it is reasonable to assume that they add in a simple 
way so that quantum mutual information is the sum of 
the quantum and the classical correlations [6, 8-10]. 

To quantify the quantumncss of the correlation con- 
tained in a bipartite quantum state Olliver and Zurek 
[3] proposed a measure for quantum correlation known 
as quantum discord and based on a distinction between 
QIT and classical information theory (CIT). A related 
quantity concerning classical correlation was proposed by 
Henderson and Vedral [8]. A recent result that almost 
all quantum states have a nonvanishing quantum discord 
[11] shows up the relevance of studying such correlation. 

At the core of the previous quantifiers of correlations is 
the one-partition (one-side) measurement on a bipartite 
system. Thus, in a general case, those quantifiers may 
be asymmetric with respect to the choice of subsystem 
to be measured. A symmetrical quantifier of classical 
correlation based on measurement over both partitions 
of a bipartite system was proposed in Ref. [12]. It was 



assumed that the classical correlation is quantified by 
the maximum classical mutual information obtained by 
local measurements over the two partitions of the sys- 
tem. An important point about these measures of corre- 
lations is related to their computation. These measures 
are based on extremization procedures over all possible 
measurements that can be performed on the subsystems 
and thus constitute a difficult problem even numerically. 
Actually, analytical solutions for the quantum discord 
were obtained recently for a certain class of highly sym- 
metrical states [13-15]. Hence, an alternative, opera- 
tional (without any extremization procedure) quantifier 
is rather desirable. 

One approach, also based on the disturbance that a 
measurement causes in the system, was used in Ref. [16] 
where several quantifiers of correlations were proposed. 
The author characterized classical states as those not dis- 
turbed by a quantum measurement process. Another in- 
teresting proposal was presented in Ref. [17]. It was 
found that, while a certain quantity related to the work 
that can be extracted from the environment using a bi- 
partite state is nonzero for all entangled states, it also 
need not vanish for separable ones, being therefore a mea- 
sure of quantum correlation. 

Besides the characterization and quantification of clas- 
sical and quantum correlations, another important prob- 
lem is the behavior of these correlations under the action 
of decoherence. This phenomenon, mainly caused by the 
injection of noise into the system and arising from its in- 
evitable interaction with the surrounding environment, 
is responsible for the loss of quantum coherence initially 
present in the system. Recently it was noted [11, 18] that, 
for a certain class of states under Markovian dynamics, 
the quantum discord only vanishes at asymptotic time, 
contrary to what happens to entanglement, which can 
disappear at finite times [19]. These results show that 
the quantumness of correlation is more resistant to the 
action of the environment than the entanglement itself. 
Although quantum discord under decoherence does not 
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exhibit sudden death, its dynamics may be very peculiar, 
exhibiting sudden changes in behavior [15]. 

Studying how decoherence affects the correlations in 
a two-qubit composite system, Maziero and coworkers 
[15] recently proposed an operational measure to quan- 
tify both classical and quantum correlations. This result 
rests on the surprising fact that, for a suitable choice of 
the noise channel, the classical correlation is not affected 
by the decoherence process, while the quantum correla- 
tion is completely destroyed. Thus, the classical correla- 
tion may be given by the quantum mutual information 
at asymptotic time [15]. 

In this article we are interested in the dynamics of 
system-reservoir correlations under decoherence. We 
consider a nonintcracting two-qubit system under the in- 
fluence of two independent environments. The most com- 
mon noise channels (amplitude damping, phase damping, 
bit flip, bit-phase flip, and phase flip) arc studied. By an- 
alytical and numerical analysis we find that, contrary to 
what is usually stated in the literature, decoherence may 
occur without entanglement between the system and the 
environment. We also show that, in some cases, the bi- 
partite quantum correlation initially present in the sys- 
tem is completely evaporated and not transferred to the 
environments, as can occur for entanglement under am- 
plitude damping as reported in Ref. [20] . 

The article is organized as follows. In Sec. II we dis- 
cuss some proposed measures of correlation: the mutual 
information in the realm of CIT and QIT, the quantum 
discord and its generalization to a "two-side" measure 
of quantum correlation. We also present a recently pro- 
posed operational measure (without any extremization 
procedure) for both classical and quantum correlations. 
A brief review of the dynamics of open quantum systems 
is presented in Sec. Ill and our results on the dynamics 
of correlations under decoherence is presented in Sec. IV. 
We summarize our conclusion and some possible avenues 
for future research in Sec. V. 



II. MEASURES OF CORRELATION 

A. Classical information theory 

In CIT the mutual information measures the correla- 
tion between two random variables A and B [21] 

I C (A:B)=U(A)+H(B)-H(A,B), (1) 

where H(X) = -X^P* lo gPz and K(A,B) = 
— J2 a bPa,b^ogp a ,b are the Shannon entropies (through- 
out this article all logarithms are taken to base 2) for the 
variable X (X = A, B) and the joint system AB, respec- 
tively. p a fi is the joint probability of the variables A and 
B assuming the values a and b, respectively, and p a = 
Y,bPa,b (Pb = Y, a Pa,b) is the marginal probability of the 
variable A {B) assuming the value a (b). 



From Bayes' rule [21], we can write the conditional 
probability 



and we also can express the joint entropy as H(A, B) = 
H{A\B)+H{B), where H(A\B) = - Y. a .bPa,b logf> a | 6 is 
the conditional entropy of the variable A given that vari- 
able B is known. Hence, the classical mutual information 
can also be expressed in terms of the conditional entropy 
as 

J C (A:B) =H(A)-H(A\B). (3) 

It is then straightforward to see that the two expressions 
(1) and (3) for the classical mutual information are equiv- 
alent [I C (A:B) - J C (A:B) = 0]. 

B. Quantum information theory 

In QIT, the extension of Eq. (1) to a bipartite quantum 
state {pab) is trivially obtained as [22-24] 

I(PA:B) = S( PA ) + S(p B )-S(p AB ), (4) 

where S(p) = — Ti(plogp) is the von Neumann entropy 
and Pa(b) — r ^ r B(A){PAB) is the reduced-density opera- 
tor of the partition A(B). It is largely accepted that the 
quantum mutual information X(pa-.b) is the information- 
theoretic measure of the total correlation (including both 
classical and quantum correlations) in a bipartite quan- 
tum state [6, 7]. 

In the context of QIT, there is no quantum extension 
of Bayes' rule (for a general state) [25]. In fact, an ana- 
log of Bayes' rule can hold only for composite quantum 
states without quantum correlation (a purely classically 
correlated system). This departure: from CIT arises from 
the nature of the measurement process in quantum me- 
chanics. Differently from the classical scenario, the con- 
ditional probability Eq. (2) in QIT depends on which 
observable is measured in the system B, since, in gen- 
eral, a quantum measurement disturbs the system. This 
implies the nonequivalence [3] of the quantum extensions 
of Eqs. (1) and (3). 

One-side measures of correlation — To obtain a quan- 
tum version of Eq. (3), let us consider a projective mea- 
surement IE- ; on the subsystem B of the composite state 
Pab- The reduced state of subsystem A after the mea- 
surement, is given by 

Pa = ^ M (1a ® nf )) pab (u ® nf >) } , 

where qj = Ttab {{lA^Uf^pAB} is the probability 
for the measurement of the jth state in subsystem B 
and 1a is the identity operator for subsystem A. For 

a complete set of projective measurements jn^ i, we 
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can define the conditional entropy of subsystem A, for a 
known subsystem B, as Sf n ( B ) j (pa\b) = J2j QjS (pa) ■ 
So, we have the following quantum extension for Eq. (3) 



J{Pa-.b) = S(p A ) - S| n ( B ) j (pa\b) ■ 



(5) 



For a quantum correlated state, Eqs. (4) and (5) are 
not equivalent. The difference 



V(pab) =Z(pa:b) - max J(pa-.b) 



(6) 



was called quantum discord by Olliver and Zurek [3] . One 
can say that Eq. (6) reveals the quantumness of the 
correlation between the partitions A and B since it shows 
the departure between QIT and CIT. We note that the 
nonclassical correlation captured by the quantum discord 
may be present even in separable states [3] . 

A quantum composite state may also have a classical 
correlation C(pab), which for bipartite quantum states 
may be quantified via the measure proposed by Hender- 
son and Vedral [8] 



C(pab) 



max 

K s) } 



S(p A ) - S 



{<>} 



(pa\b) 



(7) 



where the maximum is taken over the complete set of pro- 
jective measurements j on subsystem B. We con- 
sider only projective measurements, rather than the more 
general positive operator-valued measure (POVM) used 
in the original definition [8] of Eq. (7). In fact, Hamieh 
and cowrokers [26] showed that, for a two-qubit system, 
the projective measurement is the POVM that maximizes 
Eq. (7). For the purposes of this article we will only need 
to compute correlations between two qubits, justifying 
the restricted set of measurements. 

From the previous definitions, it follows immediately 
that V(pab) +C{pab) = I(pa-.b), as expected. For pure 
states, we have a special situation where the quantum 
discord is equal to the entropy of entanglement and also 
equal to the Henderson- Vedral classical correlation. In 
other words, V(pab) = C(pab) = 1(pa-.b)/2 [6, 8]. In 
this case, the total amount of quantum correlation is cap- 
tured by an entanglement measure. On the other hand, 
for mixed states, the entanglement is only a part of this 
nonclassical correlation [3-5]. 

It is worth mentioning that, for a general state, the 
quantum discord in Eq. (6) and also the (one-side) clas- 
sical correlation of Eq. (7) may be asymmetric with re- 
spect to the choice of system to be measured. It can be 
verified that, for states with maximally mixed marginals 
(Tta(b) Pab oc 1b(A)), D(pab) and C(pab) are symmet- 
ric under the interchange A ^ B. 

Two-side measures of correlation — Besides "one-side" 
measures of quantum (6) and classical (7) correlations, 
we can define "two-side" measures for these correlations 
[12, 27]. The classical correlation in a composite bipar- 
tite system can be expressed as the "maximum classical 



mutual information" that can be obtained by local mea- 
surements on both partitions of a composite state [12] 



K-(Pab) 



max [I c {pa:b)} 



(8) 



where I c (pa-.b) is the classical mutual information de- 
fined in Eq. (1), H(A), H(B), H(A,B) being the en- 
tropies of the probability distribution of the subsystems 
(A and B) and the composite system (AB) resulting from 

a set of local projective measurements Tl^ <g> Hj B ^ on 
both subsystems. Hence, a two-side measure of quantum 
correlation can be defined as 



Q(Pab) = Z(pa-.b) ~ IC(pab)- 



(9) 



For composite states of two qubits with maximally 
mixed marginals, we have numerically verified that the 
quantum discord (6) is identical to the two-side measure 
of quantum correlation (9) [i.e. V(pab) = Q(pab) and 
also IC(pab) = C(pab)}- 

Operational measures of correlation — Recently for a 
two-qubit system, some of us proposed an operational 
measure of quantum and classical correlations based on 
the dynamics of these correlations under decoherence 
[15]. It was shown that, under suitable conditions, the 
classical correlation is unaffected by decoherence. Such 
dynamic behavior leads to an operational measure of 
both classical and quantum correlations that can be com- 
puted without any extremization procedure. This can 
be done by sending the component parts of a composed 
state through local channels that preserve its classical 
correlation so that the quantum correlation Q(pab) will 
be given simply by the difference between the state mu- 
tual information T(pa-.b) and the completely decohered 
mutual information X [e(pa:b)] 



ab) = Z(pa-.b) - X [s(pa-.b)] 



(10) 



since the classical correlation, C(pab), present in pab is 
given by 



AB 



) =Z[e(PA:B)] 



(11) 



Here, e(pab) represents the evolved state of the system 
under suitable local decoherence channels, described as 
a completely positive trace-preserving map e(-), in the 
asymptotic time [15]. The choosing of suitable channels 
that preserve the classical correlation is the challenging 
part of this measure. Until now, this problem was solved 
only for a given class of composite states of two qubits 
with maximally mixed marginals [15]. 



III. DYNAMICS OF OPEN QUANTUM 
SYSTEMS 



Let us briefly review the theory of open quantum sys- 
tems (for a complete treatment see Ref. [28]). The time 
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evolution of a general closed quantum system is governed 
by the Liouville-von Neumann equation (we will use nat- 
ural units, such that H = 1) 

p(t) = -i[H,p(t)}, (12) 

where p and H are the density operator and the Hamil- 
tonian of the system, respectively. This equation implies 
that the evolution is unitary. However, in a realistic sce- 
nario the system of interest (5) — hereafter referred to 
as just system — interacts with its surrounding envi- 
ronment (E) (also referred to as reservoir). To account 
for this unavoidable interaction, which is often the ma- 
jor source of noise introduced into the system, we can 
rewrite the complete Hamiltonian as 

H = Hs + He + Hi, 

where Hs and He are the bare system and environ- 
ment Hamiltonians, respectively, and Hi the interaction 
Hamiltonian. While it is true that the whole system 
(S + E) still respects Eq. (12) (the density operator 
P = Pse now also includes the variables of the environ- 
ment), in general, we are only interested in obtaining an 
effective dynamic equation for the S variables. This may 
be done by taking the partial trace of Eq. (12) over the E 
variables. Then the reduced system dynamics is governed 

by 

p s (t) = -i r Tv E {[H,p SE (t)]}, (13) 

where ps — Tr# (pse) is the reduced-density operator of 
the system. This evolution is not, in general, unitary and 
leads to the phenomenon known as dccohcrencc [29]. If 
we assume that the environment is Markovian (which im- 
plies a large number of degrees of freedom) and initially 
uncorrelated with the system S (psb(O) = Ps(Q)®Pe{0), 
Pe being the reduced-density operator of the environ- 
ment), Eq. (13) can be written as a sum of operators 
acting only on the system 

p s (t) = -i[H s ,ps(t)} 

~ 7l J [Ps(t) L i L j + LiLjp s (t) 

-L jP s(t)Li]+H.c., 

where Lj is the so-called Lindblad operator and 7^ is a 
constant that depends on the specific decoherence pro- 
cess. This is the well known master equation approach 
for open quantum systems [28]. It is important to note 
that this approach depends on the perturbation theory 
for the system-environment coupling parameter, which 
implies that it is valid only in the weak coupling regime, 
i.e., when S is nearly closed. 

Although the master equation approach is widely used, 
specially in quantum optics [30], there is another way to 
treat open quantum systems, which is more appropriate 
for our purposes. We will only sketch this approach in 
what follows (a complete treatment can be found in Ref. 



[22]). The formal solution of Eq. (12) can be written in 
the form 

PSE(t) = U(t)p S E(0)lJi(t), (14) 

where U(t) is the unitary evolution operator generated by 
the total (S + E) Hamiltonian. The partial trace over the 
environment variables defines a completely positive map 
e (•) for all classically correlated system-environment ini- 
tial states [31] that describes the evolution of the system 
S under the action of the environment E 

e(ps)=Tr E {TJ(t)p S E(0)U\t)}. (15) 

The map e is a quantum operation, not necessarily uni- 
tary, mapping density operators into density operators 
and this is the reason that e is a completely positive 
map 1 . Assuming that the system and the environment 
arc initially uncorrelated [pse(0) = Ps(0) <8> p E (0)}, we 
can rewrite Eq. (15) in the so-called operator-sum rep- 
resentation 

e(p s ) = J2r kPs Tl (16) 

k 

with the Kraus operators Tk(t) = E (k\ U(t)p E \k) e act- 
ing only on the state space of system S, {|A;)_e} being an 
orthonormal basis for the environment. The Kraus op- 
erators satisfy the completeness relation ~^2 k T\.Tk = 1, 
yielding a map e that is trace-preserving 2 . This defini- 
tion of the Kraus operators is not unique. If we adopt 
a different basis to compute the trace in Eq. (15), we 
obtain a different set of equivalent operators in the sense 
that both sets generate the same dynamics for the sys- 
tem (the same operation). This can be seen from the 
fact that these two sets of operators are connected to 
each other by a unitary transformation. Moreover, it 
can be shown that, assuming a Markovian environment, 
Eq. (16) leads to the same master equation as that ob- 
tained from Eq. (13). From these considerations we can 
see that the operator-sum representation is more general 
than the master equation approach in the sense that the 
former can be applied even if the environment has only 
a few degrees of freedom. Another advantage of this tool 
is that it can be applied, in a simultaneous way, to a 
large range of physical systems, since Eq. (16) does not 
include specific details of the environment, providing us 
with a quite general dynamic equation for the system S. 

To generalize this formalism for the case in which the 
system S is composed of more than one part, we must 



1 If Aab is a positive map and pab is the density operator of the 
composite system AB, then pab = ^AB (pab) is also a valid 
density operator (all its eigenvalues are non-negative). If is 
a completely positive map, pab = I A® ^B (pab) is also a valid 
density operator. 

2 In fact, this condition can be generalized to include non-trace- 
preserving maps, such as a measurement process. In this case, 
we have Efc r I r fc < 1 [ 22 1- 
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specify which type of environment we are dealing with. 
Let us consider two types of environment (i) global and 
(ii) local. In case (i), the interaction of all parts of S 
with the same environment may lead, in principle, to an 
increase in correlations between the parts of the system 
due to "nonlocal interactions" mediated by the environ- 
ment [32] . In case (ii) , each part of S interacts with its 
local, independent environment. It is clear that corre- 
lations cannot then be increased between the parts of 
the system by interaction with the environment. For the 
case (ii), regarding N system parts and N independent 
environments, Eq. (16) immediately becomes 

e(Ps)= £ l4 , ®"< ) rf®-®r£» t . (17) 

ki,...,kN 



given that 

Use\Qs ® \0)e = 5Zr fc |Ct>s ® \k) E - 

k 

Here we will be interested only in the case of local, 
independent environments. 

IV. CORRELATION DYNAMICS UNDER 
DECOHERENCE 

In this section we will investigate the correlation dy- 
namics of a two-qubit system pab under the action of two 
local environments. The most general two-qubit state 
can be written in the form [35] 

PAB(0) = \j2 c i,M A) ®*f\ (19) 

i,j=0 

Ik) 

where a\ is the standard Pauli matrix in direction i (i = 
1,2,3) acting on the space of subsystem k (k — A,B), 
(Tq^ = 1^ being the identity operator for the partition k 
and Cjj are real coefficients that satisfy both positivity 
and normalization of pab- Our goal here is to study 
the dynamics of classical and quantum correlations as 
well as the entanglement within the possible bipartitions 
of the complete system (the system of interest plus the 
environments) under the action of several noise channels. 
We consider the most common decoherence channels (i.e., 
amplitude damping, phase damping, bit flip, bit-phase 
flip, and phase flip). 



Here 1^ is the k a th Kraus operator for the environment 
acting on subsystem a. This can be verified directly from 
the fact that the total evolution operator in Eq. (15) can 
be written in the product form U(t) = Ui(i) (g> U2(i) ® 
• • ■ <g> Ujv(i). 



The decoherence process can also be represented by a 
map in terms of the complete system-environment state. 
Let {|0)s}i w ith I — l,--- ,d, be a complete basis for 
S. Then, there are, at most, d 2 Kraus operators [33] and 
the dynamics of the complete system can be represented 
by the following map [34] 



(18) 
I 

In what follows, we will consider a system S consti- 
tuted by the two qubits A and B, each of them interact- 
ing independently with its own environment Ea and Eb, 
respectively. 

A. Amplitude-damping 

The amplitude-damping channel is a classical noise 
process describing the dissipative interaction between the 
system and the environment. There is an exchange of en- 
ergy between S and E, such that S is driven into thermal 
equilibrium with E. This channel may be modeled by 
treating E as a large collection of independent harmonic 
oscillators interacting weakly with S, as in the case of the 
spontaneous emission of an excited atom in the vacuum 
electromagnetic field (the reservoirs are at zero temper- 
ature, i.e., in the vacuum state) [28-30]. 

The action of a dissipative channel over one qubit can 
be represented by the following phcnomcnological map 
(from Eq. 18) 

|0)s|0) B -> \0)s\0)e (20a) 
|l)s|0) B -> v^|1>s|0)b + VP|0)s|l) B , (20b) 

where |0)g is the ground and |l)s the excited qubit state. 
\0)e and \1)e describe the states of the environment 
with no excitation and one excitation distributed over 
all its modes, respectively. Equation (20a) describes the 
fact that if the system and the environment start in the 
ground state, there is no dynamic evolution. Equation 
(20b) tells us that if the qubit starts in the excited state, 



ICi)s ® \o)e -> r |Ci>5 ® \o) E + ■■■ + r^K^s ® \<r - i) E 
IC2>s ® \o)e -> r |C 2 >5 ® |o) B + • • • + r d 2_ 1 |c 2 )s ® \d 2 - i) E 

\td)s ® \0)e -> r |Cd)s ® |0) B + • • • + TffiACfis ® \d 2 - 1) E 



6 



there is a probability q = 1 — p that it will remain in this 
state and a probability p that it will decay. We are using 
p to describe these probabilities as a paramctrization of 
time such that p € [0,1]. The exact dependence of p on 
time will depend, of course, on the specific model for the 
environment as well as on the system under considera- 
tion. For example, if we consider an infinite bosonic en- 
vironment interacting with a two-level fermionic system 
under the Markovian approximation, p will be a decreas- 
ing exponential function of time. On the other hand, 
if we are dealing with an (non-Markovian) environment 
with a small number of degrees of freedom, p may be an 
oscillatory function of time. The advantage of using p 
instead of an explicit function of time is the possibility of 
describing a wide range of physical systems in the same 
dynamics. 

We can get a geometrical picture of the action of this 
channel by considering the Bloch sphere representation of 
one qubit interacting with an infinite bosonic reservoir at 
zero temperature (a tractable model for the amplitude- 
damping channel). As already noted, in this case p is an 
exponentially decreasing function of time and the action 
of the channel is then to move every point of the unit 
sphere to the pole where the |0)s state is located. In 
other words, in the asymptotic limit the whole sphere is 
compressed to a single point located at the (lower energy) 
pole. 

From Eqs. (20a) and (20b), we can see that the Kraus 
operators describing the amplitude-damping channel arc 
given [22] by 



"l " 


r (fc) 


'o Vp 





> 1 1 — 






(21) 



where k = A,B labels the two distinct environments (one 
for each qubit). 

Let us start by studying the correlations between the 
various bipartitions of the whole system (S + E), assum- 
ing that the initial total density operator is given by 

PABE A E B = \{jl °^ A) \°°)eaE b , (22) 

where |00) e a e b IS the vacuum (zero temperature) state 
of the environments Ea and Eb, in which the qubits A 
and B, respectively, are immersed. The coefficients 
(co = 1), are real constants constrained in such a way 
that pabEaEb is positive and normalized. We note that 
the state of the system AB in Eq. (22) represents a 
considerable class of states including the Werner (|ci| = 
l c 2| = | C3 j = a) and Bell basis states (|ci| = |c 2 | = | C3 j = 
1)- 

Although wc can compute from Eqs. (20a), (20b) and 
(22), the total density operator, pabe a e b (p) evolved un- 
der the action of the amplitude-damping channel in a 
straightforward way, it is too cumbersome to be shown 
here. As we arc interested in the correlations be- 
tween the parts in the various bipartitions of the com- 
plete system, we will only need the corresponding re- 
duced matrices. The reduced-density operator for the 
partition AB, obtained by taking the partial trace of 
PabEaEb (p) over the degrees of freedom of the reser- 
voir pab (p) = T^e a e b [pabe a e b (p)}, in the computa- 
tional basis {|00) fc ,,|01) fcI ,|10) fcJ ,|ll) fc ,} for the parti- 
tion kl (k = A, I = B), is given by 



"(1+p) 



PAB (p) 



+ (1 




q(ci - c 2 ) 



pfc 3 

(1 -c 3 )q+(l + c 3 )pq 
q (ci + c 2 ) 




(1 





q (ci + c 2 ) 

C3 )<? + (i + c 3 ; 




pq 



q (ci - c 2 ) 



9 2 (l + c 3 )J 



(23) 



For the partitions AEa and AEg , the reduced-density 
operators read 



and 



PAE A (p) 



10 00- 

p ^/pq 

Vpv 1 





(24) 



PAE B (p) 



(l + cs)(l+pg) + l-C3 (ci-c 2 )vM 

(l-c 3 )p+(l + c 3 )p 2 

(l-c 3 )q+(l + c 3 )q 2 

{ci-c 2 )y/pq (l + c 3 )pq . 

I 



(25) 



respectively. 



Finally, for the partition EaE Bi obtained by tracing 
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out the system degrees of freedom, we obtain 



(1 



PE A E B (?) = ^ 






(ci - c 2 ) 



(ci-ca)p 

(1 - c 3 )p+ (1 + c 3 )pg (ci+c 2 )p 

(ci+c 2 )p (1 - c 3 )p + (1 + c 3 )pq 

(l + c 3 )p 2 



(26) 



Due to the inherent symmetry of the system, the den- 
sity matrix representing the partition BEb is identical 
to that for the partition AEa, thus leading to the same 
dynamics. The same symmetry is exhibited between the 
partitions AEb and BE a. 

Due to the X structure of the density matrices (23) 
through (26), there is a simple closed expression for the 
concurrence present in all bipartitions 

C(p)-2max{0,A 1 (p),A 2 (p)}, (27) 

withAi(p) = |pi4 1 - VP22P33 and A 2 (p) = \pis\- y/puPii. 
For the classical Eq. (8) and quantum Eq. (9) corre- 
lations we have no analytical expression; however, nu- 
merical analysis is possible. To this end, we will con- 
sider a Werner initial state, where c\ = c 2 = c 3 = —a 
(0 < a < 1). In Fig. 1, we show the dynamics of the cor- 
relations for the partition AB. First we note that both 
classical (JC) and quantum (Q) correlations only vanish 
in the asymptotic limit (p = 1), while the entanglement 
suffers a sudden death (SD) at a certain parametrized 
time psd [19, 36]. This can be seen directly from Eq. 
(27). On the other hand, such a system exhibits a sud- 
den birth (SB) of entanglement between the reservoirs 
(E A E B ) [20]. 

The fact that the "entanglement sudden death" be- 
tween the qubits and "entanglement sudden birth" be- 
tween the reservoirs may occur at different instants was 
first reported in Rcf. [20] and is shown in Figs. 2 and 
3. On the other hand, as we can see in Figs. 1 and 2, 
contrary to what happens to entanglement, the vanish- 
ing of the classical and quantum correlations between the 
parts of system AB is accompanied, simultaneously, by 
the creation of such correlations between the reservoirs. 
Moreover, in Fig. 1 we see that, although the entan- 
glement in partition AB disappears at a finite time, the 
quantum correlation Q vanishes only asymptotically, as 
previously noted in Ref. [18]. 

In Figs. 4 and 5 we plot the dynamics of correlations 
for the partitions AEa and AEb, respectively. From 
these figures and from the fact that pae a {p) = Pbe b (p) 
and pae b (p) = Pbe a {p), we see that each qubit has 
nonclassical correlation only with its own reservoir for 
all values of a and p. Also, in the asymptotic limit, as 
expected, all the correlations between the system and the 
reservoirs vanish due to the fact that we considered the 
reservoir initially in the vacuum state. 



FIG. 1: (Color online) Correlation dynamics for the 
amplitude-damping channel, considering partition AB for the 
Werner initial state, (a) Classical correlation [Eq. (8)], (b) 
quantum correlation [Eq. (9)], (c) mutual information [Eq. 
(4)], and (d) concurrence [Eq. (27)]. 



B. Phase-damping 

The phase-damping channel describes the loss of quan- 
tum coherence without loss of energy. It leads to decoher- 
ence without relaxation. An example of a physical system 
described by this channel is the random scattering of a 
photon in a waveguide [22]. The map that describes the 
action of this channel on a one-qubit system is given by 

|0)s|0) B -> |0)s|0) B 

|1)s|0) b Vq\±)s\0) E + Vp\1)s\1)e. (28) 

So that there is no exchange of energy between the sys- 
tem and the reservoir, only the phase relations between 
the energy eigenstates of the system are lost during time 
evolution. The Kraus operators describing the phase- 
damping channel for the qubit k (k = A, B) may be 
written as 



"1 " 


r (fe) _ 


'0 " 


,/q_ 


i 1 1 — 


s/p_ 



8 



FIG. 2: (Color online) Correlation dynamics for the 
amplitude-damping channe, considering partition EaEb for 
the Werner initial state, (a) Classical correlation [Eq. (8)], 
(b) quantum correlation [Eq. (9)], (c) mutual information 
[Eq. (4)], and (d) concurrence [Eq. (27)]. 



FIG. 4: (Color online) Correlation dynamics for the 
amplitude-damping channel, considering the partition AEa 
for the general state (22). Classical correlation (dot-dashed 
line) given by Eq. (8), quantum correlation (dotted line) given 
by Eq. (9), mutual information (solid line) given by Eq. (4), 
and concurrence (dashed line) given by Eq. (27). 



FIG. 3: (Color online) Correlation dynamics for the two-qubit 
system with the amplitude-damping channel for the Werner 
state with (a) a = 0.5 and (b) a = 0.6. Quantum corre- 
lations given by Eq. (9) for partitions AB (solid line) and 
EaEb (dashed line) and entanglement given by Eq. (27) for 
partitions AB (dotted line) and EaEb (dot-dashed line). 



Assuming the initial state (22), the evolved density op- 
erator of the partition AB, obtained by tracing out the 
degrees of freedom of the reservoirs, is given by 



FIG. 5: (Color online) Correlation dynamics for the 
amplitude-damping channel considering the partition AE B 
for the Werner initial state, (a) Classical correlation [Eq. 
(8)], (b) quantum correlation [Eq. (9)], (c) mutual informa- 
tion [Eq. (4)], and (d) entanglement [Eq. (27)]. 



£(Pab) 



1 + c 3 

1 — c 3 c + q 

c+q 1 - c 3 

. c~q 



c q 



l + c 3 . 



(29) 



where we defined c+ = C\±C2- The classical and quan- 
tum correlations present in the reduced state (29) can 
be computed analytically through the measures (11) and 



(10) and are given by [15] 



C[e( PAB )]=J2 



i + (-i) fe x 



log 2 [1 + (-l) k x\ , (30a) 



k=l 



V [e {p AB )] = 2 + X k lo §2 Afe - C Upab)] , (30b) 



k=l 
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where x — maxjg 2 |ci| ,q 2 | C2 1 , | C3 J }■ and Xk is the kth 
eigenvalue of the reduced-density matrix pab (p) [15]. 
We can verify that C [e (pab)} and V [e (pab)} are sym- 
metric under the interchange i H B. It was also nu- 
merically verified that the "one-side" measures of corre- 
lations (7) and (6) lead in this special case to the same 
values as the "two-side" measures (8) and (9), respec- 
tively. Therefore, for this state, the quantum discord 
and the Henderson- Vcdral classical correlation are suit- 
able measures of correlations. 



The correlations in the partition AEa are contained 
in the following reduced-density operator, obtained by 
taking the partial trace over the subsystems B and E B 



PAE A 



while for the partition AEb one obtains 





"1 








- 


1 














2 








1-p 


Vpv 












p . 



(31) 



P AEb (P) 



-l + q + pcs (l-c 3 )y^pq 

{l-c 3 )^pq (l-cs)p 

l + q-pc 3 {1 + C 3 )^pq 

{l + c 3 )Jpq (1 + C 3 )p . 



(32) 



The last partition we want to analyze here is EaEb, 
whose reduced-density operator of wich is given by 



PE a E„ (p) 



Aq + (1 + c 3 ) p 2 
l\fPQ 

- (1 + c 3 )m 



■~f\/pq 
ip 

(1 + c 3 )pq 
(1 + c 3 )pJpq 



(1 



-c 3 )pq 
IP 

(1 + c 3 )pJpq 



(1 



(1 
(1 



(1 



\-c 3 )pq 
c 3 )Ps/pq 

c3)pVpq 

fc 3 )p 2 . 



(33) 



r 



where 7 = 2 — (1 + c 3 )p. For these last states we must 
use the "two-side" measures of correlations (8) and (9). 
Before proceeding with the numerical analysis, let us 
look at the entanglement between the various partitions. 
Defining p^f as the partial transposition of matrix p ki 
with respect to the subsystem k [37], we can directly 
see that P t a a Ea (p) = Pae a (p), p T A % B (?) = Pae b (p), and 

Pe e (p) = PE a E b (p)- From the Peres separability crite- 
rion [37], we see that there is no entanglement between 
the subsystems A and E A m) nor between the reservoirs 
Ea and Eb, for any value of the parametrized time p. Al- 
though no bipartite entanglement was observed beyond 
that contained in the two-qubit initial state, multipartite 
entanglement is always possible. To investigate this pos- 
sibility, we consider the following system-reservoir initial 
state: 



Mi*) - -7= (\0) A |1)b - |0) b ) ® \0)e a |0> 



which is obtained from Eq. (22) by setting c\ = c 2 = 
c 3 = — 1 (Werner state with a = 1). The action of the 
phase-damping channel (28) on the above state results in 



the following asymptotic (p = 1) system-reservoir state: 

m - 77I (\°)a l 1 )^ |0> £a \1) Eb |1>a |0>b \1) Ea |0) Eb ) 

(34) 

Therefore, in the asymptotic limit, the system- 
reservoir state is a quadripartite entangled state, the 
Greenberg-Hornc-Zeilinger (GHZ) state [38]. The GHZ 
class of states is the only one that possesses irreducible 
multipartite correlations [40]. This means that the cor- 
relations in state (34) cannot be determined by looking 
at its reduced-density operators [40]. Retrieving to en- 
tanglement, there is no tripartite or bipartite entangle- 
ments in state (34) because the reduced-density operators 
are separable. This example indicates that although the 
bipartite entanglement between AEa and BEb is null, 
multipartite entanglement between all parts of the global 
system can be generated during the decoherence process. 

The next example illustrates another important fea- 
ture of the two-qubit dynamics under a phase-damping 
channel described previously. Consider a separable 
system-reservoir initial state, for example pabe a e b (0) = 
Pab(0) <8> \0) Ea (0| ® \0) Eb (0|, p AB (0) being the Werner 
state with a < 1/3. As the initial state is separable 
(and the environments arc independent) and there is in- 



10 



teraction only between the partitions AEa and BEb, no 
multipartite entanglement is generated for all values of 
p. Therefore, entanglement can only be created between 
the qubit A(B) and its reservoir Ea{Eb) due to their 
interaction. However, as shown in Fig. 6, the qubit A(B) 
does not get entangled with its reservoir for any value of 
p. Although no bipartite or multipartite entanglement 
is created during the time evolution, we do have deco- 
herence, as can be seen from the asymptotic limit of Eq. 
(29). A possible explanation for this fact is the pres- 
ence of nonclassical correlations between the qubit and 
its reservoir (see Fig. 6). Decoherence without entangle- 
ment between the system and the reservoir were noted 
earlier in the context of continuous variables [39]. On 
the other hand, when one considers a single qubit under 
phase damping, the qubit decoherence process is always 
accompanied by entanglement between the qubit and its 
reservoir. This result was verified numerically for many 
qubit initial states with nonzero coherence. 

As we can see in Fig. 7, under phase damping, the 
quantum correlation Q is asymptotically null, but the 
classical correlation JC reaches its maximum at this limit. 
Comparing this figure with Fig. 6 , where the correla- 
tions for the partition AEa are plotted, we see that the 
reduction of correlations in partition AB is accompanied 
by the creation of correlations in partition AEa- We note 
that Q between the partitions A{B) and Ea(Eb) starts 
to increase until it reaches a maximum, decreasing to 
zero thereafter and leading, in the asymptotic limit, to a 
classical correlated state between the qubits and its reser- 
voirs (the symmetry of the initial state leads to the same 
evolution for the partition BEb)- In Fig. 8 we show the 
correlation dynamics for the partition AEb where the 
same behavior as that of partition AEa is exhibited. 

From Figs. 6 to 8, we observe that the quantum cor- 
relations (including the initial entanglement) disappear 
in the asymptotic regime (p = 1), in all partitions con- 
sidered. This is not the case for the amplitude-damping 
channel, where the quantum correlations are completely 
transferred from the system AB to the reservoirs EaEb 
at p = 1. 



C. Bit-Flip, bit-phase-flip and phase-flip channels 

The effect of bit-flip, bit-phase-flip and phase- flip chan- 
nels is to destroy the information contained in the phase 
relations without an exchange of energy. The action of 
these channels on a single qubit can be described by the 
following Kraus operators 



r (fc) 
1 o 



1 
1 



4 k) = y/pal 



(fe) 



(35) 



where i = 1 (x axis) for the bit-flip, i — 2 (y axis) for the 
bit-phase-flip, i = 3 (z axis) for the phase flip (k = A,B), 
and we defined q' — 1 — p/2. 

It is helpful to get a geometrical picture by looking at 
the Bloch sphere representation of one qubit [22] . To this 



FIG. 6: (Color online) Correlation dynamics for the dephas- 
ing channel considering the partition AEa, for the general 
state (22). (a) Classical correlation [Eq. (8)], (b) quantum 
correlation [Eq. (9)] and (c) mutual information [Eq. (4)]. 



FIG. 7: (Color online) Correlation dynamics for the dephas- 
ing channel, considering the partition AEa, for the general 
state (22). Classical correlation (dashed line) given by Eq. 
(8), quantum correlation (solid line) given by Eq. (9) and 
concurrence (dotted line) given by Eq. (27). 



end, let us then consider the action of the bit-flip chan- 
nel. Owing to the symmetry of the Kraus operator (it 
is proportional to <r x ), all the points on the sphere are 
uniformly compressed on the x axis. Then states on this 
axis will be invariant under the bit-flip channel, as can 
be seen directly from Eq. (35). It is not difficult to sec 
that the actions of the other two channels are completely 
equivalent to that of the bit flip, the only difference being 
the symmetry axis. The bit-phase flip channel will leave 
invariant states on the y axis, while for the phase flip 
channel, the symmetry axis is the z axis. For this rea- 
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From these equations we directly see that P^ Ea (p) = 
Pae a (p) and p T / EB (p) = p AEB (p) implying, once more, 
from the Peres separability criteria [37] that we have de- 
coherence without entanglement between the qubits and 
the reservoirs for any parametrized time p. The last par- 
tition is given by 



FIG. 8: (Color online) Correlation dynamics for the dephasing 
channel, considering partition AEb, for the general state (22). 
(a) classical correlation [Eq. (8)], (b) quantum correlation 
[Eq. (9)] and (c) mutual information [Eq. (4)]. 



Pe a e b (p) 



(?') 




c lP q'/2 





pqi/2 
cipq'/l 






c lP q'/2 
pq'2 




Cipq'/2 



P 2 /4 



(P) = PE A E B ip)- 



(39) 
Fig- 



which also has the property P EaEe 
ure 9 shows similar behavior to that in Fig. 6 for the 
phase-damping channel. During the process of decoher- 
ence, subsystem A becomes quantum correlated (but not 
entangled) with its own reservoir and the asymptotic 
state possesses only classical correlations. 



son, wc will present here the case of the bit flip channel 
alone. For completeness, the phase-flip and bit-phase-flip 
channels are presented in the Appendix A. 

Considering the initial state (22), the evolved reduced 
density matrix for the partition AB under bit flip is given 

by 



£{Pab) 



1 



1 + c 3 q 2 













1 - 

Cl 



- c 3 g 

f c 2 q 2 




Cl 

1 - 



r C 2 <T 
c 3 <7 2 




ci - c 2 q 



l + c 3 g 2 J 
(36) 

Once more, we can use "one-side" measures of classical 
(7) and quantum (6) correlations, which can be computed 
analytically for this case [15]. They are given by Eqs. 
(30a) and (30b), but with x = max {l c i| , 9 2 l c 2| , q 2 \c%\}- 
Note that the axes y and z are continuously contracted 
by the factor q 2 while the x axes is left invariant. 

The bipartitions of subsystem A and each of the two 
reservoirs are given by 



Pae a (p) = 



q' 
o 



LVpqTz 



o 

p/2 

VpYJ? 
o 



and 



Pae b (p) 





q' 




VpYJi 
o 
o 

P/2 J 



(37) 



q' c lv /K72' 

p/2 Cl ^pYj2 

_ Cl ^pYj2 q' 
Civ^V2 p/2 . 

(38) 



FIG. 9: (Color online) Correlation dynamics for the bit flip 
channel, considering the partition AEa, for the general state 
(22). Classical correlation (dashed line) given by Eq. (8), 
quantum correlation (solid line) given by Eq. (9) and the 
concurrence (dotted line) given by Eq. (27). 



For completeness, in Figs. 10 and 11 we plot the cor- 
relation dynamics for the partitions AEb and E^Eg, 
respectively. As we can see, the dynamic behavior of the 
correlations under the bit-flip channel is essentially the 
same as that under the phase-damping channel (Figs. 6 
and 8). 



V. SUMMARY AND CONCLUDING REMARKS 

In this article we investigate the system-reservoir dy- 
namics of both classical and quantum correlations in the 
decoherence phenomenon. We considered all possible bi- 
partitions of a two-qubit system interacting with two 
local, independent environments, modeling several com- 
mon noise sources: amplitude-damping, phase-damping, 
bit-flip, bit-phase-flip, and phase-flip channels. 
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FIG. 10: (Color online) Correlation dynamics for the bit flip 
channel, considering partition AEb, for the general state (22). 
(a) Classical correlation [Eq. (8)], (b) quantum correlation 
[Eq. (9)], and (c) mutual information [Eq. (4)]. 



system exhibit some degree of correlation, including en- 
tanglement. In case (ii), the classical and quantum corre- 
lations initially present in the system are transferred over 
time to all bipartitions of the complete system, but the 
entanglement is not transferred. The bipartite entangle- 
ment contained in the system is completely "evaporated" 
by the action of the channels. While all bipartitions de- 
velop quantum correlation of separable states during evo- 
lution, at the asymptotic time all nonclassical bipartite 
correlations are null. Thus, the asymptotic state of the 
whole system (system of interest plus environment) con- 
tains only classical correlation in all bipartitions. In case 
(ii), we have decoherence without entanglement between 
the qubits and the environment, the classical and quan- 
tum correlations (of separable states) are responsible for 
the information transfer from the system to the environ- 
ment. 

Finally we note that we study here only bipartite cor- 
relations. Certainly a study of multipartite correlations 
will make a useful contribution to understanding the dy- 
namics of information in the decoherence process. An 
important future investigation will be the study of the ef- 
fects of finite-temperature environments on the dynamics 
of these correlations. Another interesting line of research 
is the dynamic behavior of the system under the action of 
a single environment, where correlations can be created 
in the system through nonlocal interactions mediated by 
the environment. 



Acknowledgments 



We are grateful for the funding from UFABC, CAPES, 
FAPESP, CNPq, and the Brazilian National Institute for 
Science and Technology of Quantum Information (INCT- 
IQ). 



FIG. 11: (Color online) Correlation dynamics for the bit flip 
channel, considering partition EaEb, for the general state 
(22). (a) classical correlation [Eq. (8)], (b) quantum correla- 
tion [Eq. (9)], and (c) mutual information [Eq. 4)]. 



We observe here two distinct behaviors for the dynam- 
ics of correlations, when the qubits are under the action 
of (i) amplitude-damping and (ii) phase-damping, bit- 
flip, bit-phase- flip and phase-flip channels. In case (i), 
all correlations (classical and quantum, including entan- 
glement) initially present in the system are completely 
transferred in the asymptotic time, to the environments. 
During time evolution, all bipartitions of the complete 



Appendix A: Correlation Dynamics for bit-phase-flip 
and phase-flip Channels 



The AB reduced-density matrix for the bit-phase-flip 
channel is given by 



Pab ip) 



' 1 + c 3 q 2 



fiiq 2 - c 2 











ciq - c 2 


" 1 c 2 l~c 3 q 2 



1 - c 3 q 2 c x q 2 + c 2 



aq 











1 + Cs<f 



the classical and quantum correlations be- 
ing given by Eqs. (30a) and (30b), with 
X = max{g 2 |ci| , |c 2 | ,q 2 |c3|}. We see that the en- 
tire Bloch sphere is shrunk onto the y axis, which is 
the symmetry axis for the bit-phase-flip channel. The 
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density operator for the partition AEa is given by and for the partition AEb, by 





r q' 










PAE A (p) = ^ 






p/2 










*vW/ 2 








P/2 



PAB B (f>) 





" 1-P/2 








-ic 2V / pq'/2 


1 





p/2 


-ic 2 ^pq'/2 





2 







q' 







ic 2y /pq'/2 








p/2 



where we directly see from the Peres partial transposition 
criterion [37] that there is no bipartite entanglement in 
these cases. The last bipartition reads 



PE a E b (p) 



(q'T 

pq'/2 c 2 pq'/2 

c 2 pq'/2 pq'2 
c 2 pq'/2 



c 2 pq'/2 



p 2 /4 



From this last equation we see once more that the reser- 
voir never gets entangled. 

For the phase-flip channel, the partition AB is 



Pab 0) = 





"1 + c 3 








c " 


1 





1-C 3 


c+q 2 





4 





c+q 2 


1-C 3 







. c~q 2 








1 + c 3 _ 



the classical and quantum correlations be- 
ing given by Eqs. (30a) and (30b), with 
X = max {q 2 \a\ , q 2 \c 2 \ , | C3 1 J- . We see that the 
symmetry axis for the bit-phase-flip channel is the z axis 
(as for dephasing) — the Bloch sphere is compressed 
onto this axis. The density operator for the partition 
AEa is given by 



Pae a (p) 



l-p/2 VSV2 
VpYJZ P/2 





and for the partition AEb is 






q' 

-VpYJz 






VpW2 
P/2 . 



P AEb (?) = 



l-p/2 c 5 ^pY]2 

c 3 ^72 p/2 

q' -cz^pT/2 

-c 3 ^/pYj2 p/2 



where we directly see from the Peres partial transposition which also does not exhibit entanglement, 
criterion [37] that there is no bipartite entanglement in 
these cases. The last bipartition reads 

' (qf c 3 pq'/2 

pq'/2 c 3 pq'/2 

c 3 pq'/2 pq'2 

c 3 pq'/2 p 2 /A 



PE a E„ (p) = 
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